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• The sensitivity of the deuteron form factor to contributions from the ex- 

cited states of the nucleon is explored using a simple model of the nucleon- 
£ — , 

, nucleon interaction which employs a tower of charged nucleon resonances. 

The model is manifestly covariant, analytically solvable, and gauge invariant. 
The consequences of this model are studied in the simplest possible frame- 
"q ■ work. We assume that all particles have spin zero and that the tower has 

only three charged members, which consist of the proton, the Roper, and a 
higher state in the vicinity of the -D13. Nucleon- nucleon S-wave phase shifts 
and the deuteron form factor are calculated using this three member tower, 
and the results are compared to similar calculations using the proton ground 
■ state only. We conclude that the deuteron form factor is insensitive to the 

presence of excited states of the proton unless those states are of sufficiently 
low mass to produce strong inelasticities in NN scattering channels. 
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I. OVERVIEW, RESULTS, AND CONCLUSIONS 



A. Introduction 



Now that first results from experiments performed at the Thomas Jefferson National 
Accelerator Facility (Jefferson Lab) are becoming available, the details of the two nucleon 
current may be explored in great detail. Some of the new experiments will probe the short 
range structure of this current, which can be interpreted as arising from the exchange of 
quarks and gluons, or alternatively, from meson exchange and the excitation of virtual 
excited states of the nucleons. In this paper we turn our attention to virtual contributions 
from these nucleon excited states. 

The effects of nucleon resonances (or excited states) are obvious above the production 
threshold; it is less obvious what effect they will have on cross sections and observables 
measured well below the production threshold. Denoting the mass of the final state in 
electron scattering by M x , and the square of the four momentum transfer by Q 2 > 0, we try 
to answer the following question: "If M x lies below the threshold for the physical excitation 
of nucleon resonances, does their existance nevertheless affect the high Q 2 behavior of form 
factors and structure functions?" To study this question we introduce the "tower of states" 
model. This is a simple, dynamical model for the study of nucleon resonances which can 
make easy, but admittedly rough, predictions of the effects that might be seen at places like 
Jefferson Lab. The goal is to use this simple model to get a qualitative feeling of where 
effects due to resonances will be large. 

Any such model must be compatible with highly successful theories of NN interactions 
[ I | - [[J which have gone before. Motivated by previous work [[| and the relativistic work of 
Tjon and collaborators [|Kj - [|L|] we use manifestly covariant dynamics. When working at 



the level of a few GeV, as at Jefferson Lab, one must use relativistic equations. The model is 
constructed so that gauge invariance is satisfied exactly. This requires the introduction of an 
interaction current which is derived from the strong kernel. Finally, the model is analytically 
solvable, allowing for simple calculations of various properties of the NN system. Many of 
these features have appeared in earlier models, but, to our knowledge, this is the first time 
the contributions from resonances and their corresponding interaction currents have been 
systematically applied to the study of the deuteron form factor using a simple model which 
is analytically solvable, covariant, and gauge invariant. 

In this first study we use the simplest possible potential (rank-one separable), with a 
minimal number of free parameters (five, at most). We fix all free parameters by fitting the 
NN phase shifts below lab energies of 350 MeV. All deuteron properties (binding energy, 
wavef unctions, form factors) are then calculated from the model. Our treatment of the 
nucleon resonances varies somewhat from earlier works. Most models have considered the 
addition of only one additional resonance, usually the A (Ref. @] is an exception, however, 
but see the discussion in Ref. WW). Such channels have been considered as add-ons to 



existing successful NN models |7j, or used to increase the applicability of such models [13 



Previous calculations of the deuteron form factor using separable models have usually ignored 
the interaction current which should accompany the calculation. The philosophy behind 
our treatment is that the resonances and their interaction currents should be consistently 
included on an equal footing from the beginning. 
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Since we wish to keep this first exploration of the contribution of resonances as simple 
as possible, we have neglected spin (all particles in the following discussions are scalar), and 
used a basic separable four-point interaction diagram as the kernel with a "tower" consisting 
of only three states. A Tabakin-style form factor, similar to the one used by Rupp and Tjon 
|I4| , is employed in order to realistically model the contributions of meson exchange without 
actually having to deal with its complexities. We compare results from a three member tower 
with those obtained from a purely elastic scattering model (a one member tower, in our 
language). Therefore, any differences between the models which include higher resonances 
and those which do not are due to the presence of the resonances. Clearly, if differences 
between these models can be seen even at this basic level, then further investigation with 
more complicated systems is warranted. 

In this section we begin with an overview of the techniques used in formulating the 
model, and then discuss the results of this work. Particular attention is paid to the effects of 
inelasticity on the tower models. At the end of this section we present our conclusions. Sec. [Tl 
includes a full development of the model, including the solution of the relativistic scattering 
equation, proof of unitarity, a proof of gauge invariance, and the development of the deuteron 
form factor. The relations used for inelastic scattering are also derived. In order to maintain 
gauge invariance, contributions from an interaction current are required, since the kernel of 
the this model is dependent on the center-of-mass four-momentum. The interaction current 
is fully derived in the Appendix, and the contributions from the interaction current prove 
to be very important to several facets of the NN system. Finally, Sec. [Ill] includes some 
specific details of the models. 



B. Background and overview 

The tower of states model is built upon the idea that, in the process of repeated scattering 
with a neutron, a proton may "sample" any one of its higher order resonances (the "tower"). 
That is, the proton and neutron interact, causing the proton to become a virtual Roper for 
a short period of time before interacting again with the neutron. In successive scatterings, 
multiple samplings of this tower of higher resonant states could occur. We explore the 
effects that such sampling will have on various properties of the NN system, particularly 
the properties of the deuteron where all inelastic contributions are virtual. 

As a brief overview of the techniques used in formulating the model, we begin with the 
most general form of the relativistic scattering equation 

M = V + / VGM, (1.1) 

Jk 

where V is the relativistic kernel, G is the propagator, and J k is specified by choosing a 
particular formalism, such as Bethe-Salpeter [IjJ, Blankenbecler-Sugar [|TB[, or Spectator 



p7| - [pCfl . The basic four-point interaction between a member of the tower and a neutron 
shown in Fig. [IJ is chosen as the kernel. The coupling strengths are denoted by the coupling 
matrix g^, where i and j denote members of the tower. We include a Tabakin-style form 
factor 
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FIG. 1. Basic tower of states interaction. This diagram forms the kernel of the spectator 
equation describing the scattering and bound states of the system. The cross on the neutron lines 
denotes that the particle is on the mass shell. 



with 
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x = x(p,P)=p ^— , (1.3) 

where p is the relative momentum between the neutron and the tower member and P is the 
center of mass four-momentum (for a detailed definition of the form factor, see Sec. HI A| ). 



As shown by Rupp and Tjon [T3J], and again here for a one-member tower called Model 1, 
the Tabakin form factor allows one to successfully model the iViV system without having to 
actually perform meson exchange calculations. Since spin is being neglected at this time, all 
particles have scalar propagators. Finally, we choose to work in the spectator formalism fl7f 



which has proven to be very useful in describing the NN system H] , |H| . For simplicity 
in this first treatment, we will treat the np system only, and discard the isospn formalism, 
so that the neutron and proton are nonidentical particles. The spectator formalism is then 
simplified if the neutral particle, i.e. the neutron, is placed consistently on its mass shell 
and excluded from the tower of states. The towers will contain the proton and its charged 
excited states only. 

These choices are not made arbitrarily. Each has been specifically chosen to reduce the 
complexity of the system in some way without losing the essential physics. By constructing 
a very basic, separable kernel, the M-matrix may be obtained analytically (see Sec. [II] for 
the full details). Neglecting spin is, of course, an automatic reduction of complexity useful 
for this first treatment of the problem. The choice of the spectator formalism motivates 
the initial inclusion of only one tower (since one particle is on shell, it is natural to restrict 
it to its ground state, although a more complete treatment can be developed later) and is 
chosen because it is easily extended to more complex cases. Finally, treating the proton and 
neutron as nonidentical particles and neglecting the distribution of charge inside the neutron 
allows us to restrict the electromagnetic interaction to the tower members; the neutron does 
not interact with the photon in the naive models considered here. 

Once the scattering equation has been solved, we obtain the S'-wave phase shifts, 5, and 
inelasticity parameters, rj = cos p, from the usual representation of the inelastic scattering 
amplitude 

-M = 5£_L i (1.4) 
2ip 1 
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FIG. 2. The basic electromagnetic coupling between tower members 

where p\ (not to be confused with the inelastic phase angle p) is the relativistic phase space 
factor given in Eq. ( |2.19| ) below. The free parameters of the models are then fit to the 



Nijmegen 1993 phase shifts below 350 MeV Excellent fits are obtained using towers 



with either one or three members. Fits using a two member tower were attempted, but 
proved to be less successful than those obtained from one- and three-member towers. The 
reason for these poor fits is under investigation; it may be related to the need for a certain 
symmetry between the coupling constants. Due to the larger number of free parameters, 
X 2 s for three-member tower models are slightly improved over those using only the proton. 
These results will be presented in Sec. |TC| . 

In order to explore the effect of the resonances in a simple NN system where they can can 
only make virtual contributions, we study the deuteron form factor. The electromagnetic 
coupling between tower members, illustrated in Fig. |2], is chosen to be proportional to 
Cij = Qje, where e is the electric charge. The full currents with electromagnetic form factors 
are constructed using the methods of Gross and Riska [[22|] (see Eq. fl2.29|) for more detail). In 
this simple first treatment we choose Qj = 5{j so the photon itself cannot excite or de-excite 
members of the tower. In order to satisfy gauge invariance an interaction current must be 
included, and the inclusion and study of this current is one of the novel features of this 
work. This current is derived from the kernel V by minimal substitution using the methods 
of Ito et. al. ||23|| . However, because the structure of the kernel is different from the one 
assumed by Ito et. al., a slightly different argument is needed, and because the derivation is 
somewhat lengthly the details are presented in the Appendix. 



C. Numerical results 

From the real and imaginary parts of M we calculate the phase shifts and perform a 
least-squares fit to fix the free parameters in the models. The resulting parameters for three 
models are shown in Table |. Model 1 consists of a tower with the proton only, Model 
3 employs a three- member tower consisting of the proton, Roper, and D 13 , and Model 31 
uses the same three-member tower, except that the mass of the Roper has been lowered to 
increase its inelastic contributions. The inelastic region has no effect on the fitting procedure, 
as only data below 350 MeV are used in the fit. 

In order to best understand the nature of the restrictions placed on these parameters, 
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TABLE I. Parameters in the tower of states models. Numbers in bold face were varied during 
the fitting procedure. The form factor parameters and all masses are in GeV; all the coupling 
constants have dimensions GeV . 



Model 1 Model 3 Model 31 



Parameter 


'So 


3 Si 


'So 




3 Si 


'So 




3 Si 


a 


0.198125 


0.155094 


0.197611 




0.157080 


0.197750 




0.150054 


P 


1.09772 


1.27872 


1.09795 




1.29825 


1.10168 




1.28994 


7 


0.171623 


0.136404 


0.172138 




0.139047 


0.172141 




0.129887 


Pc 


0.3492 


0.412689 


0.3492 




0.412689 


0.3492 




0.412689 


911 


-1516.96 


-2553.31 


-1500.89 




-2615.86 


-1504.81 




-2554.04 


912 






155.0 




406.731 


190.0 




404.563 


913 






155.0 




-105.0 


190.0 




-105.0 


922 






10.0 




-105.0 


15.0 




-105.0 


923 






272.175 




597.206 


253.928 




591.489 


933 






10.0 




-105.0 


15.0 




-105.0 


x 2 


6.97 


7.60 


5.06 




6.83 


3.66 




6.04 


mi 


0.93825 






0.93825 






0.93825 




mi 








1.44 






1.17 












1.52 






1.52 





look at the expression for tan 5. For a one- member tower (Model 1), this expression reduces 
to 

where Q\ is the integral given in Eq. (|2.7|) . In order for the phase shift to cross zero (which it 
does at around 300 MeV lab energy), there must be a zero in the numerator, and \mQ\ must 
be set equal to zero at the particular value of the center of mass energy at which the phase 
shift crosses zero. If we fix p c by this method, the zero in the numerator will be a double 
zero, and will result in the desired structure for the phase shift only if the denominator 
also has a single zero located at precisely the same center of mass energy This zero in the 
denominator of tan Si is used to determine the parameter gu. Note that the value of p c will 
be a constant during the fitting procedure (it does not depend on any of the free parameters) 
and gu will be determined actively during the fitting process. Of the five parameters which 
define Model 1, only a, (3, and 7 are are varied freely during the fit, and these are shown in 
bold face in Table |. This procedure was used for both the singlet and triplet cases. 

The careful reader will have noticed that the denominator of tan 8\ in Eq. (|1.5|) must be 
the same as the denominator of the M-matrix. Thus, this expression must also have a second 
zero in the triplet case, located at the mass of the deuteron. No attempt is made to fix the 
location of this zero, we instead determine it after the fitting is complete, obtaining the mass 
of the deuteron directly from the model. Table [Tl] shows the values of the binding energy 

TABLE II. Binding energy of the deuteron as calculated in each of the three models. The 
percent error shown is the error between the calculated value and the expected value of 2.22 MeV. 
The energies are given in MeV. 





Model 1 


Model 3 


Model 31 


Binding Energy 


2.15 


2.02 


2.47 


Percent Error 


3.2 


9.0 


11.3 
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TABLE III. Values for the phase shifts, in degrees, below 350 MeV. T LA b is in MeV. 



Nijmcgen 1993 a 



Model 1 



Model 3 



Model 31 



Tlab 


^0 


3 Si 


'So 


3 Si 




3 Si 


'So 


3 Si 


1 


62.029 


147.75 


61.559 


147.97 


60.889 


147.58 


61.2 


148.36 


5 


63.627 


118.18 


64.203 


118.27 


63.99 


117.83 


64.148 


118.54 


10 


59.956 


102.61 


60.55 


102.32 


60.498 


101.96 


60.614 


102.35 


25 


50.903 


80.629 


51.199 


80.009 


51.32 


79.786 


51.403 


79.705 


50 


40.545 


62.767 


40.675 


62.641 


40.884 


62.468 


40.961 


62.222 


75 


32.933 


51.585 


32.912 


51.852 


33.149 


51.685 


33.232 


51.435 


100 


26.783 


43.227 


26.564 


43.72 


26.812 


43.556 


26.903 


43.332 


125 


21.549 


36.46 


21.174 


37.146 


21.423 


36.985 


21.525 


36.793 


150 


16.939 


30.72 


16.371 


31.480 


16.618 


31.327 


16.73 


31.165 


175 


12.776 


25.701 


12.075 


26.536 


12.317 


26.394 


12.439 


26.258 


200 


8.943 


21.216 


8.08 


22.026 


8.3152 


21.898 


8.4471 


21.786 


225 


5.357 


17.14 


4.4249 


17.959 


4.648 


17.847 


4.7894 


17.755 


250 


1.959 


13.386 


1.031 


14.189 


1.2076 


14.095 


1.3552 


14.022 


275 


-1.3 


9.89 


-2.3423 


10.696 


-2.0353 


10.621 


-1.8673 


10.565 


300 


-4.457 


6.602 


-5.3305 


7.4381 


-5.0731 


7.3836 


-4.897 


7.3436 


325 


-7.545 


3.484 


-8.2173 


4.958 


-7.9706 


4.3246 


-7.7858 


4.3003 


350 


-10.59 


0.502 


-10.96 


1.4582 


-10.72 


1.4466 


-10.525 


1.438 



1 Reference 21] 



of the deuteron for the various models, and the percent difference from the experimental 
binding energy of the deuteron. While the these errors are large on the scale of the deuteron 
binding energy, 2.2 MeV, they result from a delicate cancellation between two large numbers 
(of the order of 2 GeV), and hence the agreement is really quite good. In light of this, all 
three models give quite satisfactory values for the binding energy. 

In the case of Models 3 and 31 there are ten parameters: six coupling constants and 
four form factor parameters. Unlike the case of Model 1, the singlet and triplet couplings 
have different restrictions placed upon them. For the singlet case, p c will be fixed using the 
location of the zero in the phase shift, as detailed above, reducing the number of parameters 
to nine. It was determined that g\\ needs to be a fitted parameter in this case, in order 
to reduce \ 2 to its optimum value. There are, however, restrictions on several of the other 
coupling constants, specifically, gi2 = gi3 and g 2 2 = #33- These four couplings are then fixed 
to constant values, reducing the number of parameters to five. 

The choice of values for certain coupling constants may seem arbitrary. It was discovered 
during initial fitting trials that the fits were most sensitive to the values of the form factor 
parameters, and rather less sensitive to most of the tower couplings. Therefore, many of the 
tower couplings were fixed at historically useful values, preserving certain symmetries which 
seemed to consistently appear in successful fits. These values, for both the form factor 
parameters and the tower couplings, are given in Table |. While the fitting procedure is 
most sensitive to the values of the form factor parameters, it also shows some sensitivity to 
the values of gn and #23- These parameters were therefore allowed to vary, while the other 
couplings were fixed such that the proportions of the couplings were maintained. For the 
triplet case, the full set of restrictions on p c and gn is once again applied; only gyi and #23 
are allowed to vary. Once again, the restriction placed on the remaining coupling constants, 
#13 = #22 = #33 , is designed to maintain the proportions of the couplings. Thus, there are 
still only five free parameters in the triplet case, but they are different from those in the 
singlet case. 

We now present numerical results for all three models discussed above, as well as our 
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analysis of the significance of these results. First are the phase shifts, shown in Table [TTT] for 
the fitting region (below 350 MeV) for the three models together with the Nijmegen 1993 
values. All three models give nearly indistinguishable fits to the Nijmegen phase shifts, with 
excellent x 2 s- Clearly, there is no significant difference between the models as regards the 
phase shifts below 350 MeV. This is as expected; substantial changes in the phase shifts are 
not anticipated until one goes above the inelastic threshold. 

After the parameters are determined, the singlet and triplet phase shifts are calculated 
out to 1.2 GeV in lab energy and are shown in Fig. |3] (the Nijmegen 1993 phases are 
also shown in the fitting region). All three models fit the phase shifts extremely well. 
More importantly, note that the phase shifts for all three models are nearly identical in 
the entire range, despite the fact that the onset of inelasticity in Model 31 is at 492 MeV. 
The inelasticities for Model 31 are shown in Fig. f|, along with the inelasticities of Arndt and 
collaborators ||24|| . We have chosen to plot p, which is defined by rj = cos p. Some readers may 
be puzzled by the shape of the singlet inelasticity. In Fig. |5] we show the inelasticities for the 
1 S'o np partial wave of van Faassen and Tjon [13j together which those of our Model 31. Ours 
is more pronounced than that of van Faassen and Tjon, but the structure is similar. Some 
deviation in both the magnitude and shape of the inelasticity is expected since the source of 
inelasticity in our model arises from the production of a zero-width nucleon resonance like 
the Roper, instead of the more realistic coupling to a continum ttNN channel associated 
with the onset of production of the A resonance. 

We have also plotted the deuteron wave function for each of the three models in Fig. || 
All the wave functions are normalized in the following manner: 

*(*) = f|. (1-6) 

Clearly, the contributions from %p 2 and ip^ are very small indeed. 

Finally, we examine the deuteron form factor, shown in Fig. [7[ Note that, in spite of the 
unrealistic simplicity of the models, all three of them describe the data |25] resonably well. 
We attribute this to the fact that each has been fitted to the 3 Si phase shifts. There is a 
small gap between the three curves most prominently seen at high Q 2 , where the average 
percent difference between Model 1 and Model 31 is 12%. This effect can only be due to the 
presence of the additional resonances in the three-member models. A small shift between 
Models 3 and 31 can also be seen. 

The separation of the curves in Fig. [7], as well as the reasonable comparison to the data, 
appears to be due entirely to the presence of the interaction current. Fig. § shows the 
relativistic impulse approximation contributions for the three models, which we denote by 
Fj. We see that they are nearly identical, except for the additional minimum just visible near 
4 GeV 2 in Model 31, and would not fit the data alone. We emphasize that only the diagonal 
terms of the interaction current have been included in this calculation of the deuteron form 
factor (recall that Qj = 5^). Inclusion of the off-diagonal terms would likely increase the 
interaction current contribution, and change the final result for the deuteron form factor. 

Clearly there is an effect due to the presence of the resonances in the three-member tower 
models, most easily seen in the calculation of the deuteron form factor, but even there it 
is rather small. However, the question must be asked, are the effects truly small, or are 
they small in this case because the inelasticities themselves are small? It is with this idea 
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in mind that, in the next section, the effects of large inelasticity on the tower model shall 
be explored. 



D. An exploration of large inelasticity 

As seen in the previous section, Model 31 has a relatively small inelasticity (but bigger 
that the actual data, as shown in Fig. yet there seems to be an effect visible in the 
deuteron form factor. We now explore this effect, specifically to discover if greater inelasticity 
will produce larger effects in such calculated quantities. A criticism frequently leveled at 
current models is that the inelasticities cannot be ignored as one proceeds to larger Q 2 . 
How large do the inelasticities need to be before one can no longer ignore their effects? It is 
known that in other partial waves the true NN system has larger inelasticities than those 
found for this model, thus it seems reasonable to explore regions of large inelasticity. 

In order to accomplish this, we choose to work backwards, in a sense. Rather than try 
to fit the known inelasticities of the NN system with a tower model, we chose a simpler 
course of action, in keeping with a first exploration of these concepts. First, a set of data 
was created, based on Model 31, with much larger inelasticity. Then, a one-member tower 
model, similar to Model 1, was fit to the low energy data of this new data set. The full set 
of calculations were then completed for both the new data and the new model. This gives us 
insight into the question of large inelasticity in the simplest and most direct way possible. 

The large inelasticity data set was created by adjusting the Roper couplings of Model 
31, that is, gu, (722, and (723, such that the inelasticity p increased to 40-50 degrees, while 
maintaining the shape of the phase shifts as much as possible. The restrictions placed on 
the couplings in Sec. fl] are also kept intact. The parameters for this set of data, called 
the LI data set, can be found in Table [TV]. This table also includes the parameters of the 
one tower model, referred to as Model 1LI, which was fit to this data set. The phase shifts 
for both the LI data set and Model 1LI are given in Table [V[ Note that the changes in the 
phase shift are very slight for the singlet case, but somewhat more pronounced in the triplet 



see Table |TTT] for the Nijmegen 1993 phase shift data). Even though the \ f° r the triplet 



fit of Model 1LI seems quite large, the fit is quite good, as we will see below. The binding 



energy of the deuteron for both the data set and the model are given in Table VI 



The inelasticities for the data set LI are shown in Fig. BL and the phase shifts are shown 



in Fig. [10]. The presence of the low mass "Roper" resonance can be seen in the kink in 
the LI data set at about 500 MeV, which is at the onset of inelasticity. Fig. |TT] shows the 
deuteron wave function for the LI data set and corresponding model. It is clear that the 
inelastic effects are much stronger than in the previous case; the ip2 and ip3 components are 
about ten times stronger than before. 

Finally, we turn to the most noticable demonstration of the effects of the resonances: 



the deuteron form factor in the large inelastic case, shown in Fig. Note the pronounced 
splitting between the model and the data, as well as the slight shift in the zero. Clearly, the 
resonances are finally having an effect. Our thoughts on this matter will be summarized in 
the next section. 
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TABLE IV. Parameters in the LI data set and Model ILL Numbers in bold face were varied 
during the fitting procedure. The form factor parameters and all masses are in GeV; all the coupling 



constants have dimensions Ge\ 



Parameter 


'So 


LI Data Set 

3 Si 


^0 


Model 1LI 


3 Si 


a 


0.197750 


0.150054 


0.197579 




0.075051 


P 


1.10168 


1.28994 


1.10175 




1.31118 


7 


0.172141 


0.129887 


0.172310 




0.058073 


Pc 


0.3492 


0.412689 


0.3492 




0.412689 


911 


-1511.44 


-1906.51 


-1534.46 




-2758.70 


912 


200.0 


1000.0 








gv-i 


200.0 


-600.0 








922 


-60.0 


-600.0 








923 


1700.0 


800.0 








933 


-60.0 


-600.0 








x 2 






0.029 




31.60 


mi 




0.93825 




0.93825 




rri2 




1.17 








m 3 




1.52 









TABLE V. Values for the phase shifts below 350 MeV for the LI data set and Model ILL 



T LAB LI Data Set Model 1LI 





'So 


3 s 1 


'So 


3 S 1 


1.0 


55.63 


140.6 


55.58 


141.2 


5.0 


61.25 


109.1 


61.26 


108.3 


10.0 


58.52 


93.78 


58.56 


92.78 


25.0 


50.06 


73.11 


50.13 


73.44 


50.0 


40.01 


57.21 


40.08 


58.82 


75.0 


32.46 


47.31 


32.53 


49.36 


100.0 


26.24 


39.84 


26.30 


41.94 


125.0 


20.94 


33.78 


20.99 


35.86 


150.0 


16.2 


28.57 


16.24 


30.51 


175.0 


11.96 


24.02 


11.99 


25.79 


200.0 


8.003 


19.89 


8.025 


21.45 


225.0 


4.384 


16.17 


4.397 


17.52 


250.0 


1.022 


12.73 


1.025 


13.86 


275.0 


-2.322 


9.559 


-2.329 


10.46 


300.0 


-5.284 


6.62 


-5.301 


7.277 


325.0 


-8.146 


3.86 


-8.174 


4.266 


350.0 


-10.87 


1.284 


-10.91 


1.428 
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TABLE VI. Binding energy of the deuteron as calculated in each of the large inelastic cases. 
The percent error shown is the error between the calculated value and the expected value 2.22 
MeV. The energies are given in MeV. 





LI Data 


Model 1LI 


Binding Energy 


2.47 


2.33 


Percent Error 


11.3 


5.0 



E. Conclusions 

The major conclusions of this paper are summarized as follows. 

1. The low energy NN data can be fit by tower models using either the proton alone 
(Model 1) or a tower consisting of the proton, the Roper, and the D 13 (Models 3 and 
31). Model 1 has three parameters in both the singlet and triplet case; Models 3 and 31 
have five parameters each. It has been noted that the fits do not depend strongly on the 
coupling constants, but rather on the three form factor parameters. Therefore, certain 
symmetries have been exploited in fixing several of the remaining coupling constants. 
In the case of Model 31, the Roper mass is set equal to 1.17 GeV in order to ensure 
that some inelasticity appears within the range of interest (T L ab = — 1.2 GeV). Only 
slight effects are seen from this inelasticity on such properties of the deuteron and its 
wave function and form factor. 

2. Three-member tower models which have small inelasticity fit the data just as well 
as models having no inelasticity. The differences between these models and a model 
with no higher order resonances are slight, and thus one might conclude that no further 
investigation is warranted. However, unless models with large inelasticity are explored, 
we cannot be sure whether the effects are small because the inelasticities themselves 
are small, or because inelasticity truly has no effect on the system whatsoever. 

3. In order to explore the effects of large inelasticity, an adjusted data set was produced 
by increasing the Roper couplings in Model 31. This large inelasticity data set was 
then fit below 350 MeV using a proton-only tower model. Much larger effects of the 
resonances are seen in this case. In this case of unrealistically large inelasticity the 
deuteron form factor is changed at large Q 2 by about a factor of 2, which can certainly 
be detected by the precision experments which will be carried out at Jefferson Lab. 

In summary, our simple model suggests that the impact of nucleon resonances on the 
deuteron form factor will be substantial only if the inelasticity in the deuteron channel is 
substantial (which is not the case in the ,S-wave channel). However, we have not explored the 
effects of spin, the coupling to the Z)-wave channel, or the excitation of nucleon resonances 
by the electromagnetic current, so our conclusions can only be very preliminary. Further 
study is certainly warranted to assess the sizes of these additional effects. 

This concludes our overview and summary of this work. The following sections will 
discuss the work in greater detail, beginning with the definition of the general tower of 
states model, and concluding with model-specific results. 
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II. GENERAL THEORY 



In this part the general framework of the tower of states model will be presented. The 
relativistic scattering equation is obtained, and an analytic solution for the M-matrix pro- 
duced. Locations of bound state poles are also discussed. The subject of unitarity is then 
discussed and relations for the phase shifts obtained. Next, the inelastic scattering matrix 
is developed and the relations for the phase shift and the inelasticity parameter derived. 
Then the gauge invariance of the model is explored. Finally, we derive the expression for 
the deuteron form factor. 

A. General description of the tower of states model 

We define the tower of states to be the proton and all of its charged excited states, 
and assume that the neutron remains unexcited by the interaction while the proton can be 
transformed into any of its excited states. The interaction which carries the ith member of 
this tower (and a neutron) into the jth member (and a neutron) is illustrated in Fig. [I] and 
takes the form 

Vij = gij fW n ,D)\f[x^ n ,D)] , (2.1) 
where x(p n , D) is defined by 

x<P»,m=£-te^-, (2-2) 

with p n the four-momentum of the neutron and D the total four-momentum. Also, gij is the 
strong coupling between the ith and jth tower members, and / is the Tabakin-style form 
factor defined in Eq. ( |1.2|) . For simplicity we will sometimes use the notation x = x(p n , D) 
and x' = x(p' n , D); when necessary the arguments will be written in full to avoid confusion. 
The variable x is unchanged by the substitution p n — > p n + \D, so that x(p n , D) = x(p, D). 
Note also that our form factor is not equal to unity when both particles are on the mass 
shell, as is the case with many form factors. 

It is important to remember that the neutron is placed consistently on the mass shell at 
all times and is not a member of the tower of states. The proton and neutron are treated as 
non-identical particles, which is possible if we limit ourselves to the np system only (which 
includes the deuteron). The coupling of the photon to the neutron will be neglected in this 
simple, begining model. The subscript o shall be used to indicate quantities associated with 
the neutron, while quantities associated with the ith tower member, such as masses, will 
be given a subscript i. However, the four-momentum of the neutron will be labeled with n 
to avoid confusion with the Oth component of the vector p, while the tower particle four- 
momentum will be written without subscripts. The quantities m\ = m = 0.93825 GeV will 
be used interchangeably. Outgoing momenta are distinguished from incoming momenta by 
the use of the prime, and internal momenta are generally denoted by k, without subscripts. 
Finally, quantities written in bold face are three-vectors. 
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B. Analytic solution for the M-matrix and bound state wave functions 



We would like to find an analytic form for the scattering matrix [E2DJ, to simplify the 
calculation of phase shifts and other such important quantities. Begin with the most general 
form for a two-body scattering equation 

M(p,p';D)=V(p,p';D) + f V(p, k; D)G(k, D)M{k,p; D) , (2.3) 

Jk 

where M is the scattering amplitude, V is the kernel, J k represents an integration over the 
unspecified internal four-momentum, and G is the two-body propagator. The exact forms 
of these terms will be discussed in Sec. |TT|. 

Assuming that the kernel is sufficiently small, perturbation theory permits us to obtain 
an iterative solution to the above equation. This solution takes the form 

M = V + J VGV + J J VGVGV + ■••+(/ Vg)"v + ■ • • , (2.4) 



also shown diagrammatically in Fig. [T3|. This may be recognized as a geometric series, which 
has a simple analytic form, 

M = (l - I VG) 1 V. (2.5) 



Note that because of the form of the kernel, each term in the sum ( |2.4|) will contain a factor 
of f(x')f(x), one form factor each for the incoming and outgoing tower member. It therefore 
becomes convenient to define the reduced M-matrix A4, such that its elements are defined 
by Mij = f(,r')fU)M, r 

In terms of matrices, Eq. (|2.5|) can be written 



M = (l-gG)- x g, (2.6) 
where the matrix Q is a diagonal matrix of propagators with elements 

Gi=\ f(k)G e , (2.7) 

Jk 



and g is the matrix of the strong coupling constants {gij} defined in Sec. |II A| and Gi is the 
propagator for the £th tower member. 

To locate the bound state poles (if any) of Eq. ( |2.5|) introduce the matrix A = (1 — gQ) 
and write the solution fl2.6p as 

where C is the cof actor matrix of A. The full matrix M is then 

M ='Mso /w ' (2 ' 9) 
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Written in this form it is clear that the locations of any bound state poles can be found 
simply by solving the equation det A = 0. 

To find the bound state wave function, recall that the integral equation for the bound 
state is 



T(p,D)= I Vip,k;D)G(k,D)T(k,D) 

Jk 



(2.10) 



where T is the bound state vertex function. This equation can also be written in matrix 
notation 



(l-VG)T = 0. 

Since the interaction is separable, the solution to this equation is 



(2.11) 



c 2 



/(*) 



(2.12) 



for as many tower members as we choose. The bound state wave function is 

<b = AfGT, 



(2.13) 



where M is a normalization constant and G is the diagonal matrix of propagators with 
elements Gg. Therefore the solution of Eq. ( |2.11| ) determines the wave function up to a 
normalization constant. 



C. Unitarity 



The most general unitarity relation, derivable from Eq. ( |2.3| ) pOfl , is 

M - M I M(G — G)M 

Jk 

= M*(Q -Q*)M , (2.14) 

where the last is the matrix version of the unitary relation. Note that the bars may be 
changed to complex conjugates since there is no Dirac structure in the problem, and hence 
the equation reduces to the following conditions on the matrix elements of M 

2i Imil% = M&G t - Gf)M tj , (2.15) 

where the implied summation over I runs over all members of the tower. 

If we are below the production threshold for any tower member except the proton, 

P, = ^^ = £>1. (2.16) 

This allows us to simplify Eq. ( |2.15| ) to 
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ImMij = p x M* x Mxy (2.17) 
For determining the phase shifts, only the expression for Mn need be considered, and 

IniMn = Pl \M n \ 2 , (2.18) 
where pi is the relativisitic phase space factor 



'W 2 - 4m2 

^ = 1 isr' (2 - 19) 

Note the similarity between Eq. ( |2.18| ) and the unitarity relation for the £th partial wave in 
nonrelativistic scattering 

Im/, = -||M 2 , (2.20) 

where K is a constant, usually equal to unity in nonrelativistic theory. Therefore, the 
standard identifications may be made, 

ft <-> -M n 

f^Pi, (2-21) 

(where the minus sign in the first relation comes from a comparison of the integral equations 
and not the unitary relation) allowing us to write 

e lSl sin 8\ 

M u = i , (2.22) 

Pi 

which implies 

cos Si sin Si 



ReMi 



Pi 

ImMn = . (2.23) 
Pi 

These can be combined to yield a single expression for Si, 

tf x = arctan f (2.24) 

This is the general expression for the elastic phase shift. In order to fit to the phase shift 
data, the real and imaginary parts of Mn must be written in terms of the free parameters. 



This will be discussed in Sec. Ill 
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D. Inelasticity 



To explore the inelastic scattering region for the tower model p6| we generalize the 
structure of the scattering matrix to include the possibility of inelasticity: 

-M U = 2— , (2.25) 

2ip l 

where rj = cos p < 1 is the inelasticity parameter. Note that elastic scattering corresponds 
to r] — 1 (or p = 0), and that as rj decreases toward zero (or p toward 7r/2) the scattering 
becomes increasingly inelastic. 

It follows from fl2.25| ) that the real and imaginary parts of M n are now 

ImM; = — (rj cos 25 1 - 1) 
Pi 

ReM/ = —n sin 2Si . (2.26) 
Pi 

Solving for the two unknowns, rj and Si, gives 

piImMn + 1 



V 



cos 2(5! 
ReMj 

IniMn + 



ReMn 

tan25i = . (2.27) 



Pi 

We will return to these expressions for the phase shift and inelasticity parameter in Sec. 



E. Electromagnetic current and gauge invariance 

In this section we describe the one and two body currents, and show that the model is 
exactly gauge invariant. 

Charge conservation requires that the coupling between the photon and a given tower 
member be of the form = Qje, where e is the electromagnetic charge of the proton and 
Qj is a scaling factor which must equal unity if z = j. The off-diagonal terms are not 
constrained by charge conservation, but the simplest way to insure that the model is gauge 
invarinat is to require them to be be transverse, that is, q^j^ = 0. Finally, we require that 
the diagonal component of the current satisfy a Ward-Takahashi identity [^jj, P^j , which 
may be expressed in a compact fashion, 

= ~ e (p' 2 ) ~ G J' Ml 5 v » ( 2 - 28 ) 

where p and p' are the initial and final momentum of the tower member, respectively. This 
restriction is satisfied by the following form of the electromagnetic current of the tower 
member: 

%tf,p) = e^(Q 2 ) (p» - + e l3 S l3 {^-A , (2-29) 
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where P = p' + p and the form factors Fij(Q 2 ) represent the structure of the current due to 
the true (nonpointlike) nature of the tower members. To avoid kinematic singularities we 
require 

F«(0) = 1 

F ij (0) = 0. (2.30) 

Except for these requirements, the form factors need not be further specified at this time. 

Because the strong form factors / depend on the total four-momentum, we expect them 
to generate an interaction current. Using a technique similar to that of Ito, Buck and Gross 
||23|| , and assuming that Qj = Sij for simplicity, leads to the following interaction current 
(see the Appendix for details) 

where 



-ft p n -(D + D')(D 2 + D' 2 )Y (2.32) 
and Z fti is obtained from Z^ by replacing p n by p f n . Contracting into this current gives 

INT ( 1 

%3% (p'mPn, D>, D) = e|Vy(K,^; D') - V^pn, D)| . (2.33) 



Recalling that the bound state is generated by the infinte sum given in Eq. ( |2.4| ), it can 
be shown that the electrodisintegration process is described by the four diagrams shown in 
Fig. [L4]. We will show that these diagrams, together with the conditions ( |2.28|) and ( |2.33|) , 
are sufficient to ensure that the electrodisintegration process is gauge invariant. 

Note first that all four diagrams have the same incoming state, namely 



iB 



(p n ,D) = G e (p)Te(p n ,D), (2.34) 



where <pf is the bound state wave function of a neutron and the £th member of the tower 
of states. The final states, however, are different. For diagrams 0(a) and |H|(c) there is no 
final state interaction and the final state is a plane wave state with 

$ jpw (p' n ,Pn;D') = 2E8 ij {2ixf5\p' n -p n ) , (2.35) 



where E = ym^ + p^. Diagrams 0(b) and |i~4|(d) include the final state interactions. In the 
notation of the figures, this is 

4scatt(?4Pn;£>') = Mijip'^P^D^Gjip'). (2.36) 
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These two processes can be combined by introducing the total final state wave function 

<Hj(Pn> Pn, D ') = Mj pw(Pn> Pn] D') + (pfj scatt {p' n , p n ] D') 

= Mijip^Pn, D^Gjip') + 2E8 ij (2tt) 3 5 3 (K - p n ) . 



Using this notation, the four diagrams shown in Fig. 14 can be written 



{J")i = j j0^^P'n^D')j^p\p)<Pf{p n ,D) 

+ / ^)HEE» ?M,Pn; D') jf T (p,p"; D', D) 0f (pfi, D) , (2.37) 

where summation over repeated indicies is implied. 

Now, to demonstrate that the current ( |2.37| ) is gauge invariant we contract it with 
and use the relations ( 2.28|) and ( |2.33| ) to obtain 



= -e / T^^tij(P>n;D') {GJ 1 (p') - Gj\p)}<pf(p n ,D) 



(2tt) 3 2£ 

f d 3 Pn d Z p" n - f , 
+ e / ~, Tfi ViiXPmPn, E) ) 

x\v je (p n ,p';;D') -V jt (jp n ,&D)}<j>f(j&,D) . (2.38) 



We can reduce the expressions Gj l {p)4>f and <f>{jGj 1 {p') which occur in the first term of 
Eq. ( |2.38D by using the fact that 4>L and 4>f satisfy the spectator equation: 



y ^3 



G-\p)ct>f{p^D) = J J&jp ViJfrJj D)<j> B m {jpl D) , (2.39) 



and similiarly for the final state. Using these equations to rewrite the first term of Eq. ( |2.38| ) 
yields 



%(J% = ~e I " 6 J^, {timiPn, Pn, D')V mj (p^ Pn , D>) <f)f ( Pn , D) 

-UjiP'mPn, D')V jm (p„,K; D) ^(pl D) 



r d 3 p n d?p" n - f , 

+ e / Z Tfi <PiAPniPn\ D ) 



d 3 p n d?p" n - f u 
(2nfAEE" 

x {V jt (pny n \ D') - V jt ( Pn ,p' n] D^fipl D) 
= 0, (2.40) 

where the cancellation follows from interchanging the p n and p" n integrations in the first term 
and relabeling the sums. This completes our proof that the current is gauge invariant. 
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F. Deuteron form factor 



In this section we discuss the calculation of the deuteron current given by the two dia- 
grams shown in Fig. The deuteron form factor, Gq, is defined by the relation 

J£ = eG c {Q 2 )(D + D'y . (2.41) 

The form factor describes the structure of the deuteron; the current for a structureless point 
particle is obtained from this equation by replacing Gc — * 1. 
The current for the processes shown in Fig. [T5] is 

Jd= f 7Z^ff(P»>WM>P)<t»j(Pn,D) 



(2vr) 3 2£ 

+ / ^)HEE> % D '^r^Pn, D\ D)4>f(j> n , D) , (2.42) 

where the one-body and interaction currents are the same as those given in Eqs. (]2.29|) 
and ( |2.31| ). When the polarization vectors e M of the photon (including the longitudinal 
polarization vector associated with a virtual photon) are contracted into these currents, 
any terms containing q 11 will vanish. Using the notation shown in the relativistic impulse 
approximation (RIA) diagram given in Fig. [15| leads to the replacement = (D+D')' 1 —2p%, 
and the effective one-body current is simply 



eF(Q 2 ){(D + Dy-2p^6 lJ . (2.43) 



Since the interaction current does not contain any factors of g M , it remains as given in 
Eq. (HO). 

We now extract an explicit formula for the deuteron form factor from Eq. ( j2.42j) . First, 
note that the integrand of the integral 

h=J 4>f (Pn, D>) <P?( Pn , D) , (2.44) 

is an even function of the three vector p„ in the Breit frame, and hence in this frame 



K = I ifivn, D') K. ■t'Uvn, D) 



Pn -(D + D') 
(D + D>) 2 



(p n ,D). (2.45) 



Since both sides of this equation are covariant, this expression holds in any frame. 

The interaction current term can also be simplified in the same way. Since the initial 
and final states are identical, the integrand of the interaction integral is symmetric under 
the interchange p n <-> — p' n in the Breit frame [which transforms x(p' n , D') into x(p n , D) and 
x(p' n ,D) into x(p n ,D')\, except for the terms in Z 11 and Z' 11 which are linear in the three- 
vector parts of p^ and p'g. These terms are odd under this interchange, and thus vanish, 
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leaving only the time components of p£ and p'£- The interaction current can therefore be 
written as the product of two integrals 



INT — % e 9ij 



<f>f( P ' n ,D')f[ X (p' n ,D')] 



'(Pn,D)Z 



, . f[x(p n ,D')]-f[x( Pn ,D)]- 



(2tt) 3 2£^ v ™' ' \ x(p n ,D')-x(p n ,D) J 
where, for the deuteron form factor with D 2 = D' 2 = Mj, Z^ simplifies to 
Z» = (D + D'Y 1 
= (D + D'Y 



(2.46) 



(Pn ■ D) 2 + (p n ■ D') 2 \p n .(D + D')} 



2Mj 



M 2 d (D> + D) 2 



(2.47) 



Each of the two integrals in (|2.46| ) is covariant, and therefore can be evaluated in any frame. 
The first of these two integrals depends on D' only, and is therefore a constant (which is 
obvious when it is evaluated in the rest frame of the outgoing deuteron). Denoting this 
constant by Cj 



d 3 P ' n 



(p' n ,D')f[x(p' n ,D')} 



(2tt) 3 2£' 

and using (|2.45 ) and (|2.46| ) gives the following result for the deuteron form factor 



(2.48) 



G C (Q 2 ) = F(Q 2 

-\-2gij Ci 



d 3 p r , 



(27r) 3 2E 



\Pn,D>) 



2 Pn -(D + D') ] B 
(D + D') 2 \ 
-ff[x( Pn ,D')}-f[x( Pn ,D)}\ 



x(p n ,D') - x(p n ,D) 



(2.49) 



Note that the only the first term (the RIA) depends on the nucleon form factor, F(Q 2 ). 
The second, exchange term, is independent of this factor and depends only on the strong 
form factor /. As we will discuss in the next section, this term actually goes to a constant 
as Q 2 — ► oo. 

This concludes the discussion of the general theory of the tower of states. In the next 
section we will apply this to the specific models under consideration. 



III. TOWER OF STATES MODELS 

In this section the detailed form of the equations appropriate for a three-member tower 
are presented. Then the details of the equations needed to calculate the phase shifts and 
inelasticity parameter for a three-member tower above the inelastic threshold of the second 
tower member are presented. In the final subsection the behavior of the interaction current 
contribution to the deuteron form factor is more carefully examined. 



20 



A. Specific structure for a three-member tower 



Using the spectator formalism, the generic integrals in Sec. [11 B| become 

d 3 k 



(2vr) 3 2£ 



(3.1) 



where E = \jvn\ + k 2 (where k 2 is the square of the three- moment a in these formulae, not 
to be confused with the four- momenta), and the elements of the propagator matrix ( |2.7| ) in 
the rest frame of the NN system, where D = (W, 0), are 



d 3 k f 2 [x{k,D)] 



(2%) 3 2E E 2 -(W - E) 2 ' 



(3.2) 



The momenta are labeled in Fig. [Tj| Here E { = y mf + k 2 and, in the rest frame, 

{a 2 + k 2 ){p 2 c - k 2 ) 



f[x(k,D)} 



( 7 2 + k 2 )(p 2 + k 2 y 



f(k 2 )- 



(3.3) 



If W < m + rrii (where m = m n as discussed above) , the integrand of Eq. ( |3~2| ) is real 
for all k and the integral is real. When W > m + mj the integrand has a pole in the range 
of integration located at k = k l p , where 



\ 



'W 2 — m 2 + 
2W 



— mj 



(3.4) 



This pole gives the integral an imaginary part. The threshold for the onset of this complex 
behavior, in terms of the kinetic energy of the nucleon in the iViV laboratory system, is 



(m Q + mj) 2 — 4m 2 



LAB 



2m n 



(3.5) 



These thresholds are given in Table [VI]] for each of the channels discussed in Sec. |. 



TABLE VII. Thresholds for the onset of 



inelasticity for the various channels. 



channel mass 


LAB 


in GeV 


in MeV 


0.93825 





1.17 


492 


1.44 


1138 


1.52 


1344 
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For energies in the fitting region (Tlab < 350 MeV) only Q\ will have an imaginary part. 
(The case where more than one of the integrals has an imaginary part will be discussed in 
Sec. |U1B| .) For all energies Q\ can be written 

Gl = V l — W(2E-W) +tlmGl > (3 ' 6) 
where V denotes the Cauchy Principle Value of the integral, and the imaginary part is 



_ 4 m 2 

Imgi = vr V — ° f(kl p ). (3.7) 

The real part can be evaluated numerically. 

Now return to Eq. ( |2.9| ) and evaluate the NN scattering amplitude, which in the tower 
notation is M\\. At energies below the production threshold for either of the resonances, the 
imaginary part of M\\ can come only from Qi, which contributes only to the denominator 
det(l — gQ) = det A, and recalling that the cofactor matrix of A is denoted by C, we obtain 

M "- Re(de t y/t(de t V ( * ,)/W - ^ 

Rationalizing the denominator gives the real and imaginary parts of M n , and hence the 
phase shift in the elastic scattering region depends only on the det A 

For a three member tower, the real and imaginary parts of det A are 

Im(det A) = - \g xx - (g n g 33 - gf 3 )G s - (fl-iito - flf 2 )^2 + B g Q 2 G?!\ Im(/i 
Re(det A) = 1 - g 33 g 3 - 5-22^2 + (#22#33 - gl^GzGz 

- S-ii - (jguQaa ~ 9n)Gs ~ (fl'iito - $2)81 + B (J Q 2 Q^ Re(?i , (3.10) 

where B g is 

B g = (911922933 ~ 9ll9 2 23 - 9n933 + ^913912923 - 9X3922) ■ (3.11) 

The deuteron vertex function for a three-member tower, 

/ 

(3.12) 

is obtained from Eq. Q2.11|) . The solution for c 2 and C3 in terms of c% (which must be 
determined by normalization) is 




= £12(1 -#33^3) + 913923 g c 

(1 - 022^2) (1 - #33^3 ) - 023^2 " ' 



C = fl^ 1 - 022^2) + 912923 g c , g 

3 (1 - 022^2) (1 - 033^3) - 023^2 
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where the integrals Qi must be evaluated at W = Mo- Normalizing the wave function to 
unity at k = 0, the general expression for the deuteron wave function for a three- member 
tower is 

= = C - Mp(2m °7 Mg) (3.14) 

*f(0) Cl (ro? - ml + 2M D ,Jml + B- Ml) <* 2 vl 

Finally, we evaluate the deuteron form factor in the case of the three-member tower. The 
integrations in Eq. (|2.49 ) are easily carried out if the momenta in Fig. [15] are evaluated in 
the Breit frame, where 



p = D — p n 
P=D'-p' n 

g = (0,q) (3.15) 
D= (D ,-|q) 



D' = (Do, iq 



with D = yj + q 2 /4. We use the standard notation Q 2 = —q 2 = q 2 . For the nucleon 
form factor, F(Q 2 ), we use the dipole approximation, 



i ^ 2 



F (Q ) = (t^TotTq 5 / ' < 316 » 

where Q 2 is in GeV 2 . 

In the next subsection we extend the results for M n to the inelastic region. 



B. Inelasticities in the three-member tower 

Since one of the tower models (Model 31) has a significant inelastic region, in this sub- 
section we develop the relations needed to calculate the phase shifts and inelasticities above 
the production threshold for particle 2. 



In Section [III A| , the determination of the real and imaginary parts of M was considerably 
simplified by realizing that only terms containing Q\ could contribute to the imaginary 
part. In the inelastic region Q2 will also have an imaginary part. Breaking the numerator 
and the denominator of Eq. fl3.8| ) into their real and imaginary parts and rationalizing the 
denominator gives 



_ Im(C T ^) 11 Re(det A) - Re(^) n Im(det A) 
' ' 1 [Re(detA)] 2 + [Im(detA)] 2 ' ' 



ReM n = Re^) 11 Re(detA) + Im^) 11 Im(detA) 

[Re(detA)] 2 + [Im(detA)] 2 V ' K ! V ; 



where 
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Re(C T g) u = g n {l - g 33 G 3 ) + g 2 13 Q 3 - DReQ 2 
lm(C T g) 11 = -DlmQ 2 
Re(det A) = 1 - 5-33^3 - g u ReGi - g 22 RaQ 2 

+(911933 - gl 3 )GzRzGi + (922933 - s^^Re^ 

+ (911922 - #? 2 )(Re£iRe£ 2 - Im£ilm£ 2 ) 

— det g ^(Re^iRe^ — lmQilmQ 2 ) 

Im(det A) = -g 11 lmQ 1 - g 22 lmQ 2 + (g u g 33 - gl 3 )Q 3 lmQ 1 + (g 22 g 33 - gl 3 )Q 3 lm.g 2 
+ (911922 - gl 2 )(ReQ 1 lmQ 2 + Im£iRe£ 2 ) 

- det g Q 3 (Reg 1 lmg 2 + Im^Re^) , (3.18) 



and 



D = 



and 



#11022(1 - g^Gz) + gwg\ 3 Qz - #12(1 - g^G?) - 2g 12 g 13 g 23 g 3 + g{ 3 g2iGz 
det g = gng 22 g 33 - gng 23 - gf 2 g 33 + 2g l2 g l3 g 23 - g\ 3 g 22 . 

C. The interaction current contribution to the deuteron form factor 

We have already derived the specific formula for the three-member tower deuteron form 
factor in Section [III A| . The analytic form of the interaction current contribution, specifically 
in the limit as Q 2 approaches infinity, merits closer examination. It shall be shown that this 
limit reveals a rather unusual property of the interaction current contribution: it is constant 
as Q 2 becomes very large. 

Recall from Eq. (|2.49|) that the interaction current contribution to the deuteron form 



factor is 

r ( n 2 \-9n r f d3pn ^(n r>\? ( f ^ Pn > ~ f ^ Pn > D ^ \ 

<hc{Q ) - 2 9ij a j ^-^0, (p n , d)z y x(pniDt) _ x{pniD) ) 

= 2g ij C i l j {Q 2 ), (3.19) 



where Z and Ci are as defined in Eqs. ( g!47p and (|2~48"D . Using the fact that (j)f(p n ,D 



NGj(p)f [x(p n ,D)] [recall Eq. (|2.13| )1, Eq. ( |3.19|) gives the following expression for 1j{Q 2 ) 

This integral can be evaluated in any frame (it is covariant) and it is convenient to evaluate it 
in the incoming deuteron rest frame. In order to express the functions in Ij in the incoming 
deuteron rest frame, all of the momenta must be boosted from the Breit frame to this new 
frame. This boosting can be accomplished by the use of the matrix A, which carries D to 
rest: 



A 



/ D \ l /D M / D \ (M D 



\ 1 n I \ i_ n / \ l n I \ n 

\- 2 < l / \ 2? v o J \- 2 Q/ \ u 




(3.21) 
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where q is the magnitude of q (and will be taken to be in the z direction). Using this boost 
A, the outgoing deuteron momentum becomes 




Ml 



h 2 



M 



D 



(3.22) 



qD 



and x(p n , D) = —p 2 n which is clearly independent of q. Therefore the wave function <fif(p n , D) 
is also independent of q. However, 



x(p n , D') 



m. 



Ml 



E[M D + 



DoPnQZ 



2M, 



D 



M, 



D 



q 1 



4Ml 



(E — p n z) 2 as q — > oo 



(3.23) 



where z = cos 8, with 6 the angle between p n and z. Hence 

/ [x( Pn , D')} - f [x( Pn , D)} 4Mlf [x( Pn , D)} 
x{p n , D') — x(p n , D) q\E-p n zf 

because f(x) — > as x — > oo. Similarity, Z becomes 



(3.24) 



z E ' 1 



1 



2Ml 



2Ml 



E l 



Q 



2Ml 



D 

— PnQZ 



M\ 



D 



\2D\E - D p n qzf 
AM%Dl 



8Ml 



(E — p n z) 2 as q — > oo . 



(3.25) 



Thus the large ^dependencies of Eqs. fl3.24Q and ( |3. 25| ) cancel, and the limit of 1 as q 
approaches a constant 



oo 



UQ 2 



Af 



d 3 Pr, 



(2tt) 3 2£ 

d 3 p n 



Gj{p)f [x{p n , D)] 



\M% 



(E - p n z) 



AMlf[x( Pn ,D)} 
q\E-p n zf 



2M 2 D J (2n) 3 2E 



G 3 {p)f 2 [x(p n ,D)) 



2M 2 D Gj ' 



^3.26) 



where Qj is evaluated at D 2 = Ml. This means that the interaction current part of the 
form factor also goes to a constant. Calculating this constant from Eqs (|2.48 ) and ( 3.26 ) 
gives a value of —1.12 x 10~ 4 , which agrees very well with the numerical results obtained for 
the interaction current contribution at very large Q, as shown in Fig. pj]. We believe that 
this rather unusual feature of the interacton current comes from the fact that we are using 
a point-like four-point interaction as the kernel. At large Q the separable form factor does 
not prevent the basic point-like structure we have chosen from being seen. 
This concludes the discussion of the specifics of the tower models. 
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APPENDIX: DERIVATION OF THE INTERACTION CURRENT 

In this appendix we derive the interaction current given in Eq. (|2.31|). While the method 



is similar to that used in Ref. p3|| , there are certain differences which will be detailed here. 

The interaction current will be obtained by minimal substitution, which is most easily 
done in coordinate space. The interaction is chosen to depend on D and p n ; the coordinates 
conjugate to these momenta will be labeled r and r d . From 

p ■ r + p n ■ r n = D ■ r + p n ■ r d , (Al) 

it follows that r d = r n — r. Using these conjugate coordinates, the strong interaction in 
coordinate space can be written as (factors of g are neglected for the moment) 

r rl 4 r>' rl 4 r> rl 4 D' rl 4 D 
A = / a Pn a P" a Jf d U e i(D'.r'+ P > n .r' d ) e -i(D.r +Pn .r d ) f^')f( x ) S \D' - D) 

J (27TJ 12 



d A p' n d 4 p n d A D 



(2tt) 12 

&(y,r' d ,r d ), (A2) 



where y = r'—r, and x' and a; were defined in Eq. ( |2.2| ) (with x' obtained from x by replacing 
p n with p' n ). Note that in coordinate space, the strong interaction is not separable, as it is 
in momentum space. 

Particle 2 is neutral, so dependences of the potential on p n and p' n 

pi = -i-A- ft = -i-A- (A3) 
9r dtl dr' dlx 

do not generate any current. Particle 1 is charged, and the coordinate associated with it is 
y. Therefore, from minimal substitution the dependences on generate 

D» = -i- > -i— + e e"e iq - y . (A4) 

oyy. oy^ 

Applying this to Eq. ( |A^ ) generates an interaction current proportional to 

A(y; r' d , r d ) — > A(y; r' d , r d ) +ee" T^y; r' d , r d ) , (A5) 

where e e M T M is the first order (in the charge e) change to A which arises from the substitution 
). To find T, this change must be evaluated from 
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5A = ee»T,(y;r' d ,r d ) 



D 



f 



a 



x\p' n ,-i-^- + ee^e^' 



f 



X I Pni 1 



_d_ 



+ ee»e iq - y 



- f 



4- -4 



X I Vni 1 r% 



(A6) 



where 



/„-/ 



d 4 Dd 4 p' n d 4 p 



(27T) 



12 



To evaluate this, assume that / can be expanded in a Taylor series 

f( x ) =Y / c n x n , 

n 

define Z^ by 



(A7) 



(A8) 



xip n ,-i— + efe™ 



5 



(A9) 



and work out the difference x n — x n arising from the minimal substitution. Since 8 X — X X 
does not commute with x, this difference must first be written 



(5x)x m p n ,-i 



(A10) 



Now take into account the fact that each appearance of d/dy will act on all factors of y to 
the right of it. This includes the factor of exp i(D ■ y + p' n ■ — p n ■ r^) in the definition of 
SA where this expansion of the form factor will eventually be used. Therefore, factors of x 
occurring to the right of the factor S x will depend on D only, while factors of x occurring 
to the left of the 5 x term will depend on D + q = D' once the expansion is inserted into 
Eq. (|A6|). Hence the change in x n under minimal substitution to first order in e is 



E 



„n— 1— m 



X 



(p„, D + q) 



ee^Z^e tqy 



x m (p n ,D) 



= ee»Z, e iq y £ x™" 1 "" 1 ^, D + q)x m (p n , D) . 

m 

Using the algebraic identity 



(All) 



n-l 
m=0 



n—l—m^m 



a n -b n 
a — b 



(A12) 



Eq. ( [A 11| ) becomes 



~ r n_ r n_ pf » 7 j q -y X n (p n ,D + q)-X n (p n ,D) 

x(p n ,.D + g) -x(p n ,D) 



(A13) 
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Therefore, the change in / under minimal substitution is 

f[x(p n ,D + q)]-f[x(p n ,D)} 



6f = ee»Z^e^y 



x(p n ,D + q) -x(p n ,D, 



(A14) 



Note that Eq. ( Aq) has the form /'/ — /'/, where the tilde indicates the use of minimal 
substitution. To first order in e this can be written 



f'f - f'f = f'f - f'f + f'f - f'f 
= Sf'f + f'8f 
= 5f'f + f'5f. 



(A15) 



Thus, Eq. 



becomes 



5 A = J d A D '- y+p " r *- pn - rd \ |e Z^ 



, f[x(p' n ,D')}-f[x(p' n ,D)} 



x(j/ n ,D')-x(j/ n ,D) 



f[x(Pn,D)] 



+ee»Z fl f[x(p' n ,D')} 



^f[x(p n ,D')]-f[x( Pn ,D)] 



x(p n ,D') -x(p n ,D) 



(A16) 



Recalling that 5 A = ee^T^, we can unfold the integral to give 



■INT 
J n 



e9 \ Z » x(p> n ,D')-x(p> n ,D) t^ D )\ 



-Z,f[x(p' n ,D')} 



.J[x(p n ,D')]-f[x(p n ,D)} 



x(p n ,D') - x(p n ,D) 



(A17) 



which is the one-member tower version of Eq. ( p.31| ). 

To complete the calculation we must find Z^. We begin from its definition in Eq. (|A11|) 
and construct it directly. Using the definition of x, Eq. ( |2.2| ), assuming that x may be 
expanded in a power series in e, using the notation 



= -i 



d_ 



—%- 



d_ 



(A18) 



and keeping only those terms which are first order in e, the change in x due to minimal 
substitution is 



x — x 



(p n ■ Df + (p n ■ Df 



eD 2 
1 



eD 2 



[iPn-D) (Pn-5) + ( Pn -5) ( Pn -D)] 
+ {A (Pn ■ Df , 



1 -{6,D} + {6,D} ' 



4£> 4 



where {A, B} = AB + BA, and 



(A19) 
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{5, D} = {e • D, e iqy ) = e ■ D e iqy + e iq ' y e • D 



(A20) 



is the summetrized expansion of D 2 — D 2 . We have expanded these functions symmetri- 
cally, using anticommutators of non-commuting operators, since every factor of D contains 
a derivative which may not necessarily commute with all of its surrounding terms. The 
parameter A is introduced to provides a continuous choice between two different ways of 
symmetrizing the expansion of the D~ 2 term; it is not clear a priori which of this con- 
tinuous range of choices is correct. (If the operators all commuted, the result would be 
independent of the choice of A.) Extracting the common factor of & e iq ' y [remembering that 
factors of D contain derivatives which act all the way to the right when they are substituted 
into Eq. ([A 61)1, gives 



x — x 



Pn, \pn-(D + D')\ 



1 1 

(D + D'), \(p n ■ D) 2 + (p n ■ D') 2 

1 - A 



x 



1 1 



+ 



2D 2 D> 



(A21) 



Therefore, 



( 1 



\2D' 2 

+ (D + DY\(p n -D) 2 + (p r , 



2D 2 
D') 2 



1 1 



+ 



1 - A 
2D 2 D' 2 



(A22) 



All that remains is to determine A. To do this, look at Eq. ( |A22| ) in the case of elastic 
scattering (D 2 — D' — M|,). Contracting with q gives 



( Pn ■ D'y - (p n ■ Dy 



M 2 D 



x(p n ,D') - x(p n ,D) . 



(A23) 



This is the result needed for the total current to be conserved [this result leads directly to 
Eq. ( |2.33|) 1. Hence, current conservation requires that this should be the result even when 
the scattering is not elastic. In that case, 



(Pn ■ D') 2 
+ {Pn ■ D) 2 



1 



2D' 
1 



2D 5 



2-h 
2 

2 



2D 2 
1 



2D 



ft 



2 A,+ 
?A1 + 



XD' 2 


XD 2 




4D' A 


XD' 2 


XD 2 ' 


AD 4 


4D> 4 



(A24) 



Clearly the choice of A = will reduce the expression to the desired form, reduces to 
Eq. (|2.32|) , and the interaction current (|A17| ) is completely specified. The generalization of 
the interaction current to the case with more than one tower member is straightforward; the 
result is given in Eq. (|2.31| ). 
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FIG. 3. The singlet and triplet phase shifts to 1.2 GeV in 
lab energy for Model 1 (solid line), Model 3 (dashed line), and 
Model 31 (dotted line). The Nijmegen 1993 phase shift data 
are shown up to 350 MeV (Ref. 21) (black circles). Note that 
the phase shifts in both cases are nearly indistinguishable. 
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FIG. 4. Inelasticities for Model 31. The shape of the 1 «S'o 
inelasticity, while more pronounced than others, is not un- 
usual. The black dots are the inleasticities of Arndt and 
collaborators (Ref. 24). 
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FIG. 5. The 1 S'o inelasticities of van Faassen and Tjon 
(Ref. 12) (solid curve) and the Model 31 (dashed curve). Note 
that the shape of the Model 31 inelasticity is similar to that 
of van Faassen and Tjon, though more pronounced. 




33 




-1 .5 - " 

_2 o r 

0.1 0.2 0.3 

k (GeV) 

FIG. 6. Deuteron wave functions for Model 1 (solid 
curves), Model 3 (dashed curves), and Model 31 (dotted 
curves). In the lower graph, the bold face curves are the 
ip2 components, while the normal curves are the ip3 compo- 
nents. These wave functions are dimensionless due to the 
manner in which they are normalized. 
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FIG. 7. The deuteron form factor for the three models. 
The curves are as in Fig. 6. The black circles are data from 
Ref. 25. 
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FIG. 8. The relativistic impulse approximation contribu- 
tion to the deuteron form factor. The curves are as in Fig. 7. 
Note the presence of an additional minimum for Model 31. 
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FIG. 9. The singlet and triplet inelasticities for the LI 
data set. 
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FIG. 10. The phase shifts to 1.2 GeV for the Model 1LI 
(solid line) and the LI data set (black circles). Note that 
the kink in the triplet phase shift for the LI data set occurs 
exactly at the inelastic threshold. 



37 



0.1 0.2 0.3 

0.2 

0.1 

0.0 
-0.1 
x -0.2 
-0.3 
-0.4 
-0.5 
-0.6 

0.1 0.2 0.3 

k (GeV) 

FIG. 11. Deuteron wave functions for Model 1LI (solid 
curves) and LI data (broken curves). In the lower graph, the 
dashed line is the curve for -i/^o, while the dotted line is tp20- 
These wave functions are dimensionless due to the manner in 
which they are normalized. 
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FIG. 12. The deuteron form factor for the large inelastic model and data set. The curves are 
as in Fig. 11. 
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FIG. 13. Diagrammatic form of Eq. (2.4). 
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FIG. 14. The four diagrams required for the calculation of the electrodisintegration process. 
The large boxes with photons coupling to them are interaction currents, the oval with a V is a 
bound state vertex, and the circle with an M is a scattering matrix. 
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FIG. 15. The diagrams needed to calculate the deuteron form factor. Diagram (a) is the 
relativistic impulse approximation and diagram (b) includes the interaction current. 
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FIG. 16. The definition of the momenta for the bubble diagram. It is these momenta which 
must be used in the integrals Q%. 
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FIG. 17. The interaction current contribution to the 
deuteron form factor in the case of Model 1, plotted out to 
very large values of Q. The dashed line shows the calculated 
value of the constant which the interaction current should 
approach as Q becomes large. 
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